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Ionization of the Earth’s atmosphere by sunlight forms a complex, multi¬ 
layered plasma environment within the Earth’s magnetosphere, the inner¬ 
most layers being the ionosphere and plasmasphere. The plasmasphere is be¬ 
lieved to be embedded with cylindrical density structures (ducts) aligned along 
the Earth’s magnetic held, but direct evidence for these remains scarce. Here 
we report the hrst direct wide-angle observation of an extensive array of held- 
aligned ducts bridging the upper ionosphere and inner plasmasphere, using 
a novel ground-based imaging technique. We establish their heights and mo¬ 
tions by feature-tracking and parallax analysis. The structures are strikingly 
organized, appearing as regularly-spaced, alternating tubes of overdensities 
and underdensities strongly aligned with the Earth’s magnetic held. These 
hndings represent the hrst direct visual evidence for the existence of such struc¬ 
tures. 


Corresponding author: S. T. Loi, Sydney Institute for Astronomy, School of Physics, The 
University of Sydney, NSW 2006, Australia, (sloi5113@uni.sydney.edu.au) 
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1. Introduction 

The plasmasphere is a toroidal region within the Earth’s magnetosphere that is hlled 
with plasma from the ionosphere, which is a photoionized layer near the surface of the 
Earth [Goldstein, 2006]. The existence of held-ahgned ducts in the plasmasphere is widely 
accepted as the explanation for ground-based detections of whistler-mode waves, which 
may be guided along such ducts [Ohta et ai, 1996; Bakharev et ai, 2010]. Wave-particle 
interactions that accelerate and precipitate particles into the atmosphere occur efficiently 
within these ducts, so their existence aids the removal of energetic particles from the 
magnetosphere and thereby influences global magnetospheric-ionospheric coupling and 
energetics [Sonwalkar, 2006]. Duct properties have been studied mainly through indirect, 
ground-based observations of whistlers from both natural sources (lightning) and man¬ 
made VLF transmitters [Singh et ai, 1998; Ganguli et al, 2000; Glilverd et al, 2003]. Since 
the guiding structures cannot be detected independently of the whistlers, observations 
that use lightning-generated whistlers are necessarily intermittent, and those that use 
man-made whistlers are spatially restricted. Direct detections of plasmaspheric ducts 
have been achieved by satellite in-situ measurements [Sonwalkar et al, 1994], but these 
and also indirect satellite-based measurements such as sounding [Darrouzet et ai, 2009; 
Yearby et al, 2011; Karpachev et al, 2014] are heavily restricted by spacecraft orbital 
patterns. Duct morphologies have thus been difficult to establish. 

The Murchison Widefield Array [MWA; Lonsdale et al, 2009; Tingay et al, 2013] is a 
low-frequency (80-300 MHz) radio telescope located in Western Australia at 26°42'12"S, 
116°40'15"E (geographic), and a geomagnetic latitude of 38.6°S (Mcllwain L parameter 
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[Mcllwain, 1961] of 1.6). Although chiefly intended for studying the cosmos, it is also an 
exquisitely sensitive probe of the near-Earth plasma. Novel design concepts include its 
wide field-of-view (FoV) and high-cadence imaging capability. A multitude of unresolved 
celestial radio sources (mostly radio galaxies and quasars) back-illuminate the plasma, 
allowing one to probe fluctuations at high spatial completeness by measuring the angular 
distribution of source refractive shifts. Since interferometers measure angular positions 
using baseline phase differences, they are insensitive to the constant offset component 
of the electron column density (“total electron content”, or TEC), which adds an equal 
phase to all antennas and cancels out. Rather, the refraction-induced angular shift A6 
of a radio source depends on the transverse gradient V_lTEC in the TEC towards the 
source, and is given by [Thompson et al, 2001] 


Ae 


1 1 

-^VxTEC 

QTT^ eoUle V 


( 1 ) 


Here e and are the electron charge and mass, Cq the vacuum permittivity, and v the 
radio observing frequency. The negative sign indicates that the direction of refraction is 
towards decreasing TEC. 


2. Observations and Results 

We analyzed MWA data comprising 46 snapshots, each integrated for 2 min at a fre¬ 
quency of 183 MHz over a 30.72 MHz bandwidth, of a 30° x 50° (N-S x E-W) patch of sky 
near zenith. These were recorded over a 1.5-hr interval on 2013-10-15 between 1346-1517 
UTC (pre-midnight local time), in which the telescope was tracking celestial coordinates 
RA = 0°, Dec = —27° (J2000). We measured the refractive shifts of the ~1000 celestial 
sources visible in each ~30°-wide snapshot by fitting to the positions of the same sources 
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seen in each snapshot and computing offsets from the time-averaged positions. In more 
detail, the source-finding software Aegean [Hancock et al, 2012] was first used to lo¬ 
cate intensity peaks (candidate radio sources) in the images. We subjected the candidate 
sources to a number of quality restrictions and then cross-matched the remainder with 
the NRAO VLA Sky Survey (NVSS) catalog [Condon et al, 1998], a published database 
of celestial radio sources, retaining only those candidates with counterparts in NVSS. 
Cross-matching was performed by identifying the nearest NVSS source to the candidate 
source within the search radius (3 arcmin) and associating it to the latter. 

We took the reference position of a source, uniquely identified by its NVSS catalog name, 
to be the average of its celestial coordinates (measured by Aegean) over all snapshots. 
Then for each snapshot, we computed the angular offset vector of each source from its 
reference position. The distributions of these offsets (equivalently the TEC gradient field) 
are shown in Fig. 1 for several snapshots (see Movie SI for an animation of the full 
dataset). Strikingly organized bands of arrows are visible, with a characteristic band 
separation of ~2°. These reflect a spatially oscillating gradient field, implying alternating 
bands of underdensities and overdensities. The bands are tightly aligned with the Earth’s 
magnetic field. Taking the divergence of the vector field reveals the duct structure more 
prominently. This is shown in Fig. 2 for a representative snapshot (see Movie S2 for an 
animation of the full dataset). 

The measured offsets display a proportionality (Fig. 3a) consistent with Equation 
(1), where A is the observing wavelength. Source motions exceed measurement errors 
and are often larger than the width of the telescope point-spread function (119arcsec at 
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183 MHz). Quasi-sinusoidal oscillations of radio sources are apparent; the time series for 
one source, with errors, is shown in Fig. 3b. Note that celestial sources drift across the sky 
much faster than do the structures, which are almost stationary above the MWA. Their 
characteristic TEC gradient (a measure of the amplitude), obtained from the mean angular 
offset using Equation (1), increases over time and roughly doubles in 1.5 hr (Fig. 3c). This 
is comparable to the growth timescales predicted and reported for whistler ducts in the 
literature [Singh et al, 1998; McCormick, 2002]. 

We estimated the altitude, motion and orientation of the structures by using the MWA 
as a stereo camera. The MWA consists of 128 receiving elements (“tiles”) spread over 
an area about 3 km across (Fig. 4a). By dividing the array east-west and imaging using 
the two groups of tiles separately, we measured a parallax shift (Fig. 4b) that implies 
a characteristic altitude of 570 ± 40 km (L = 1.8), averaged over the interval and the 
FoV. This is within the L-range previously inferred for similar structures [Jacobson and 
Erickson, 1993; Hoogeveen and Jacobson, 1997]. 

A parallax analysis is only meaningful for structures that are localized in altitude. How¬ 
ever, the MWA measures the spatial derivative of a column-integrated quantity, and so is 
sensitive to irregularities over a wide range of altitudes. An argument against the struc¬ 
tures being extended in altitude, and thus for the applicability of the parallax analysis, 
comes from the observation that the bands are still prominent up to zenith angles of 
C ~ 20°. If they were extended in altitude (like vertical sheets) then oblique lines of sight 
would pierce multiple sheets, causing the features to blend into one another and wash 
out. Clear bands in the V_lTEC vector held exist out to C ~ 20°, implying that the 
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irregularities cannot be extended in altitude by more than cot C ~ 3 times the inter-band 
separation, i.e. that their horizontal and vertical cross-sectional widths are comparable. 
This further supports the cylindrical duct interpretation, as opposed to vertical sheet- 
like structures. The contribution to the TEC from the smooth plasma above and below 
the irregularity layer cancels away when the horizontal spatial derivative is taken, and 
so the situation is observationally equivalent to a thin layer of plasma embedded with 
irregularities that is bounded above and below by vacuum. This supports the validity 
of the parallax technique, and explains why we obtained a well-dehned altitude upon 
measurement. 

Computing the parallax separately over the hrst and second halves of the interval es¬ 
tablishes a downwards drift from 720 ± 90 km to 470 ± 40 km with time. Splitting the 
data spatially into northern and southern halves of the FoV reveals larger altitudes to the 
north. This is consistent with the steep magnetic inclination at the MWA site (—60°), 
and implies that the structures stretch between ~400-1000km within the FoV, likely ex¬ 
tending above this. Thus they bridge the topside ionosphere and inner plasmasphere, and 
so the transition from a neutral- to plasma-dominated medium. These properties are con¬ 
sistent with the altitudes that whistler ducts are expected to occupy at night [Sonwalkar, 
2006]. The physical spacings of the tubes are then 10-50 km, and associated TEC fluc¬ 
tuations between about 0.1 and 0.7TECU (ITECU = 10^®elm“^), again in agreement 
with inferred properties of whistler ducts [Sonwalkar, 2006]. 
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3. Discussion 

3.1. Previous Capabilities 

While the hypothesized existence of held-aligned density ducts accounts for a large body 
of experimental data [Ohta et al, 1996; Singh et ai, 1998; Carpenter and Smith, 2001], it 
has not been previously possible to directly verify that large, cylindrical density structures 
greatly extended along the Earth’s magnetic held lines exist, or to track their motions 
precisely and continuously over regional scales. Satellite in-situ observations [Sonwalkar 
et al, 1994] can only measure densities at single points, and cannot instantaneously probe 
the regional-scale density distribution. Topside sounding data and whistler spectrograms 
have been interpreted in terms of signal guidance by held-aligned ducts [Darrouzet et al, 
2009; Yearly et al, 2011], but such methods rely on assumptions about wave propagation 
[Singh et al, 1998], and detections are only possible under suitable conditions. They also 
suher from spatial restrictions imposed by satellite orbital paths and the placements of 
ground-based receivers. With a model for the Earth’s magnetic held, radial motions of 
ducts and their L-shell values can be deduced from whistler spectrograms, but zonal drifts 
cannot be similarly established [Saxton and Smith, 1989; Carpenter and Smith, 2001]. 

Claims of regularly-spaced, regional-scale, held-aligned density structures with prop¬ 
erties similar to those seen here have been made based on observations by the Very 
Large Array (VLA) radio interferometer [Jacobson and Erickson, 1993; Helmboldt et al, 
2012a, b] and the Los Alamos plasmaspheric drift radio interferometer [Jacobson et al, 
1996; Hoogeveen and Jacobson, 1997]. Although these studies could measure zonal drifts, 
identihcation of the structures relied on htting plane waves to sparsely-sampled measure¬ 
ments, a model whose suitability could not be independently verihed. The assumption 
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that the irregularities were held-aligned, together with a model for the Earth’s magnetic 
held, were required to infer their altitudes. Many of these assumptions are unnecessary 
here. The sampling completeness of the MWA surpasses that of the VLA and Los Alamos 
instruments by a factor of ~100, making it the hrst radio telescope capable of imaging the 
plasma. This has allowed us to visualize for the hrst time the regional-scale, held-aligned 
nature of these density structures, and to demonstrate a clear and remarkable spatial pe¬ 
riodicity in their layout, without any prior assumptions about their morphology. Unlike 
these and other ground-based studies, our parallax technique to establish the altitudes 
requires no knowledge of the Earth’s magnetic held. By removing the reliance on wave 
transmission for detection of the ducts, they may be monitored passively and continu¬ 
ously for long durations, including periods when structural distortions or imperfections 
may prevent them from sustaining wave propagation. 

3.2. Formation Mechanism 

Plasmaspheric density irregularities appear more often during geomagnetically dis¬ 
turbed periods, but the underlying formation mechanisms are not well established 
[Hoogeveen and Jacobson, 1997; Darrouzet et ai, 2009]. Geomagnetic conditions dur¬ 
ing the interval were mildly unsettled with a Kp index of 2, the observations taking place 
in the recovery phase of a moderate storm (maximum ring current of Dgt = —45 nT near 
04 UTC, 24-hr maximum Kp of 4). Theoretical studies indicate that ripples and undula¬ 
tions in the plasmaspheric electron density may be produced by an interchange instability 
that preferentially occurs just after periods of strong magnetospheric convection [Sazykin 
et al, 2004; Buzulukova et ai, 2008]. Such ripples have been observed on a global scale in 
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satellite observations [Goldstein, 2005]. Although below the spatial scale of the grids used 
in those simulation studies, smaller duct-like structures such as those seen in our data may 
form if this mechanism continues to operate on smaller spatial scales. We note that field- 
aligned structures of lesser prominence and multiplicity often appear in MWA data under 
quieter geomagnetic conditions. They are visible on around half of all nights inspected so 
far, most of which (with the exception of the current dataset) were not obtained during 
periods of storm-time activity. 

Organized azimuthal structures in the outer plasmasphere have been attributed to E x B 
convection under the influence of ultra-low frequency (ULF) standing waves [Adrian et al, 
2004]. Although the spatial scales of ULF waves greatly exceed those observed here, small- 
scale density structures can grow by flux-tube interchange under the imposition of a much 
larger-scale electric field [Rodger et al, 1998]. The steady downward transport established 
above is consistent with E x B convection for a westward electric held of ~lmVm“^, a 
reasonable value for periods of substorm activity [Carpenter and Smith, 2001]. Quasi- 
periodic, held-aligned structures called plasma bubbles form through a related process, 
where travelling ionospheric disturbances (TIDs) trigger the Rayleigh-Taylor instability 
[Buhari et al., 2014]. However, plasma bubble formation is only viable near the equator 
where near-horizontal geomagnetic held lines can support the plasma against gravity, 
and not at mid-latitudes. The huctuations are unlikely to be TIDs themselves, which 
have a dispersion relation [Shiokawa et al, 2013] grossly inconsistent with the spatial and 
temporal properties measured here. Furthermore, TID-like huctuations routinely appear 
in MWA data, and these exhibit markedly diherent properties from those seen here. 


DRAFT 


April 27, 2015, 12:28am 


DRAFT 



X - 12 


LOT ET AL.: REAL-TIME IMAGING OF DENSITY DUGTS 


Spatially-varying electric fields in the plasmasphere may cause localized compressions, 
squeezing plasma down selected flux tubes [Park and Banks, 1974], Secondary ionization 
from particle precipitation may also produce flux tube-selective densifications. Although 
natural precipitation events are largely conhned to the auroral zones, the powerful VLF 
navigation beacon at North West Cape (~400km north of the MWA) is known to cause 
electron precipitation concentrated towards the south [Parrot et al, 2007]. A Fourier 
analysis of our data reveals the existence of structures drifting in opposite directions. 
Components with larger spatial frequencies are observed to drift eastward, while those 
with smaller spatial frequencies drift westward. We measure drift speeds within 2% of 
corotation, consistent with these being in the inner plasmasphere where corotation is 
relatively strict [Galvan et al, 2010]. The regular spacing of the structures could have 
arisen from a Kelvin-Helmholtz instability driven by shear in the plasma, and thereafter 
amplihed through an interchange instability. 

4. Conclusion and Outlook 

Despite the seeming ubiquity of held-ahgned ducts in magnetospheric systems (e.g. lo- 
Jupiter [Imai et al, 1992] and the Sun [Duncan, 1979]), self-organization processes in 
plasmas have been difficult to isolate observationally. Our results demonstrate that wide- 
held radio telescopes such as the MWA are powerful quantitative tools for studying their 
formation, dynamics and morphology. Radio telescopes differ fundamentally from many 
approaches for probing the ionosphere and plasmasphere in that they measure density 
gradients rather than absolute density. Their insensitivity to the constant offset compo¬ 
nent of the TEC makes them particularly suitable for studying plasma density distur- 
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bances/irregularities and permits them to probe regions above the peak electron density 
with ease, since they are not signihcantly shielded by the denser underlying plasma. 

The ground-based feature tracking and altitude triangulation capabilities we have 
demonstrated here offer valuable opportunities for the real-time, regional-scale monitor¬ 
ing of inner magnetospheric structures and dynamics on a near-continuous basis, uncon¬ 
strained by the limitations of spacecraft orbits or the propagation of whistlers. Unlike 
ground-based whistler observations, which only allow for radial motions to be measured, 
the MWA can track both radial and horizontal motions, thereby allowing bulk plasma 
drifts to be characterized in three dimensions. Its potential to perform 3D reconstruction 
of density structures can provide empirical constraints on the plasma distribution both 
along and across magnetic flux tubes. This may be useful, for example, for specifying 
appropriate boundary conditions in simulation studies of global plasma flows [Tu et al, 
2006]. 
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Figure 1. The vector fields of celestial source offsets for the first (a), middle (b) and last 
(c) snapshots, overplotted with the geomagnetic field lines (black solid lines) from a single L- 
shell. Field line positions were computed based on the Australian Geomagnetic Reference Field 
model [Geoscience Australia, 2010] at an altitude of 600 km. Assumed to be locally straight and 
parallel, their apparent convergence arises from a perspective distortion due to their steep (—60°) 
inclination. The center of each arrow marks the time-averaged position (used as a reference point) 
of a celestial radio source, and the arrow represents its displacement vector. Arrow lengths have 
been scaled to 50 times the actual displacement distance. Blue and red denote arrows with 
positive and negative x components (westward and eastward), respectively. North is up and east 
is to the left as per the astronomical convention, and x = y = 0° marks the location of the zenith. 
Units along the x and y axes correspond to angular distance across the MWA FoV (1° -H- 10 km 
at 600km altitude). 


Figure 2. The vector divergence plot for the middle snapshot (i.e. the divergence of the 
vector field shown in Fig. lb), with the duct structure clearly visible. Positive divergence regions 
(relative overdensities) are represented by red-yellow hues, while negative divergence regions 
(relative underdensities) are represented by blue-cyan hues. White represents either values close 
to zero (within the FoV) or NaN (outside the FoV). Geomagnetic field lines have been overplotted 
as black solid lines, as for Fig. 1. 
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Figure 3. (a) Wavelength (A) dependence of the size of the offset for five bright radio sources, 

whose flux densities are shown in the hgure legend (1 Jy = 10“^® W Hz“^). Each set of points 
represents measurements taken in a single snapshot of time. Straight lines represent best hts to 
measurements at four different frequencies within the 30-MHz observing bandwidth, constrained 
to pass through the origin, (b) East-west component of the displacement of one of the sources in 

(a) as a function of time. Different colors represent measurements at four independent frequencies. 
Dashed lines indicate the width of the telescope point-spread function. Error bars in (a) and 

(b) represent position fitting errors, (c) Time dependence of the average magnitude of the TEC 
gradient over all sources in a snapshot, for images formed using the full 30-MHz band. 


Figure 4. (a) Layout of the 128 MWA tiles, with hollow square and circle markers highlighting 

the two groups of tiles used for measuring the parallax shift. Central-core tiles (crosses) were 
excluded to maximize the effective baseline (~890m), leaving 37 tiles for the east half and 35 
for the west, (b) Histograms of the parallax shifts measured for the east and west halves, each 
from 92 independent measurements (46 snapshots x 2 instrumental polarizations). The relative 
displacement of the two histograms represents a parallax of 5.3 ± O.darcmin. 


DRAFT 


April 27, 2015, 12:28am 


DRAFT 



0 

5 

0 

5 

0 






-15 -10 -5 0 5 10 


xn 






UTC 2013-10-15 15:17:59 



-5 


0 


5 





15 

10 

5 

0 

-5 

-10 

-15 -10 -5 0 5 10 


UTC 2013-10-15 14:31:11 





E-W offset (arcmin) 



2.2 


2.4 


2.6 


3 


3.2 










E-W offset (arcmin) 















Counts 






North (km) 



East (km) 



